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We consider the possibility of constraining the lifetime of radiatively decaying dark matter par- 
ticles in clusters of galaxies inspired from generic axions of the Kaluza-Klein type. Such axions 
have been invoked as a possible explanation for the coronal X-ray emission from the Sun. These 
axions, or similar particles, can be produced inside stars and some of them remain confined by the 
deep gravitational potential of clusters of galaxies. Specifically, we consider regions within merging 
galaxy clusters (Abell 520 and the "Bullet Cluster"), where gravitational lensing observations have 
identified massive, but baryon poor, structures. From an analysis of X-ray observations of these 
mass concentrations, and the expected photon spectrum of decaying solar KK-axions, we derive 
lower limits to the lifetime of such axions of r « 10 23 sec. However, if KK-axions constitute less 
than a few percent of the dark matter mass, this lifetime constraint is similar to that derived from 
solar KK-axions. 

PACS numbers: 14.80.Mz, 95.35.+d, 95.85.Nv, 98.65.-r 



I. INTRODUCTION 

The particle nature of the dark matter remains a mys- 
tery. Whereas cosmological observations have deter- 
mined that there is approximately six times as much 
dark matter as there is baryonic matter, we still have 
not identified which particle or particles constitute the 
dark matter. 

An apparently unrelated issue is "the solar corona 
problem" : In order to explain the quiet solar X-ray spec- 
trum, one may need to invoke new physics, such as intro- 
ducing massive axions of the Kaluza-Klein (KK) type 
P, 0. The standard axions with the remaining rest 
mass window in the sub-eV range [3| live far too long 
to solve this problem. KK-axions appear in the frame- 
work of large extra dimensions [3, Q but they are not 
considered necessarily to be dark matter candidates. In 
this work a generic example of relatively massive radia- 
tively decaying axions has been used. In the large ex- 
tra dimensions scenario, it is assumed that only grav- 
ity propagates in the higher-dimensional space, while 
the standard- model fields are confined to our (3+1)- 
dimensional subspace. Since axions are singlets under the 
standard-model gauge group, they could also propagate 
in the higher-dimensional space. As a result of compacti- 
fication, the higher-dimensional axion field is decomposed 
into a KK tower of states with the mass spacing of order 
1 / R, where R is the compactification radius of the extra 
dimensions ||. All KK excitations have the same cou- 
pling strength to matter. A source of axions will emit all 
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KK states up to the kinematic limit. Another interesting 
feature of the higher-dimensional theory is that the ax- 
ion mass m a , corresponding to the lowest KK state, may 
decouple from the PQ scale and may be determined by 
the compactification radius, m a = 1/(2R). 

We will here address the possibility that part of the 
dark matter in galaxy clusters are like KK-axions which 
are produced in normal stars and accumulated near their 
birth place over the lifetimes of these stars [H, H, 0] • Part 
of the so produced axions are confined by the gravita- 
tional potential of the host galaxy clusters, where they 
stay until they eventually decay spontaneously. These 
KK-axions are predicted to decay into two photons, pro- 
ducing a broad characteristic spectrum which happens to 
peak just where current astrophysical X-ray instruments 
are most sensitive, i.e. in the few keV range. 

We will consider X-ray observations of merging clus- 
ters, from which weak gravitational lensing observations 
have provided strong evidence for a matter component 
which is non-baryonic. Assuming that parts of this dark 
matter is of the KK-axion type, we can derive the first 
cosmological constraints on the lifetime of this kind of 
dark matter particles. 

II. RADIATIVELY DECAYING PARTICLES 

All of the massive KK-axions produced in the Sun as 
a possible explanation for the quiet solar X-ray spectrum 
1], can decay radiatively, and from the assumed produc- 
tion mechanism, we can calculate the resulting photon 
spectrum. For solar KK-axions [l[ the peak energy lies 
at approximately 4.5 keV (red line in Figure []}, and for 
E~ t > 2.0 keV the derived spectrum is well represented 
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FIG. 1: The reconstructed quiet Sun X-ray spectrum (black, 
upper curve) [8| and the expected flux from decaying massive 
axions of the Kaluza- Klein type (red, lower curve) 

by the following expression: 

F der {E^ > 2.0 keV) = 2.2 x io 8 £- 8 - 1 e -( 32 -°/' B -?) (1) 

The shape and the energy of the maximum of the spec- 
trum are tightly bound by solar physics processes. The 
normalization of the spectrum depends on the coupling of 
the KK-axions and can essentially be treated as a free pa- 
rameter. Assuming the Sun to be a "typical" star in cos- 
mos, which is a reasonable assumption, Equation [T] can 
be applied to KK-axions produced in stars and confined 
by the gravitational dark matter well of clusters of galax- 
ies. Solar X-ray observations made with the YOHKOH 
mission fit the general feature of the radial distribution 
of the massive solar axion model, except at low photon 
energies (below 4 keV) , where the bulk of the quiet and 
not-quiet solar X-ray spectrum is emitted (see black line 
of Figure 1, and discussion in ref. [10]). 

III. OBSERVED CLUSTERS OF GALAXIES 

Observations with the Chandra X-ray Observatory of 
the two clusters of galaxies Abell 520 and the "Bullet 
Cluster" (also known as 1E0657-558) were analyzed. Ba- 
sic data on the clusters and the observations are specified 
in Table 1 for a cosmology with Ho = 0.71 km/s/Mpc, 
fijif = 0.3, Q\ = 0.7 [Til] - The data were retrieved from 
the public Chandra archive and processed using standard 
data reduction methods with CIAO version 3.3 @. 

As seen from the X-ray images in Figure [5] both clus- 
ters are merging systems and the Bullet Cluster is show- 
ing a prominent bow shock feature [l2l . fl3| ] . In such 
disturbed systems the hot X-ray emitting intraclustcr 
plasma (which is the dominating baryonic component 
in clusters) is displaced from the stars and the dark 
matter, leaving mass concentrations practically devoid of 
baryons. These so-called "dark matter blobs" are ideal 
environments for dark matter studies [HI, 0, EH since 




FIG. 2: The X-ray image of Abell 520 (above) and the Bullet 
Cluster (below) overlaid the gravitational potential from weak 
gravitational lensing (red, Bullet Cluster contours from 
and the dark matter blob region (white), and the reference 
region (green). 



they have a high dark matter density and a low contam- 
ination of X-ray "noise" from baryonic matter. 



IV. DATA ANALYSIS 

For both clusters, two regions are selected for the anal- 
ysis. One region is covering a dark matter blob with high 
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mass but low X-ray emission (white in Figure [2]) . The 
other region is a reference region with the same size and 
form and of similar X-ray flux, but much smaller mass 
(green in Figure [5]). For comparison to the results of 
the region covering the dark matter blob of the Bullet 
Cluster is chosen identical to their SUB region i.e. as the 
large circle centered on the mass peak with the smaller 
circle centered on the shock front subtractred. For both 
clusters, there is a mass contrast of an order of magnitude 
between the dark matter blob region and the reference re- 
gion in which the baryons are the dominating source of 
X-ray emission. 

The dark matter in the Milky Way halo also con- 
tributes to the signal |l6[, and we must first remove 
such a contamination. Hence, before extracting the spec- 
tra, a background region with small mass and no sig- 
nificant X-ray emission, is chosen from the same obser- 
vation and subtracted. In this way any possible signal 
from the galactic halo is subtracted, simplifying the ex- 
pected signal from the decaying dark matter as only one 
source redshift has to be considered. However, the X- 
ray flux from the Milky Way (which is proportional to 
the ratio of the Milky Way halo mass within the con- 
sidered fields of view divided by the mean halo distance 
squared) is negligible compared to the cluster X-ray flux, 
i.e. M halo /D 2 halo « M blo b/D 2 blob . 

A convenient "basis model" consisting of a thermal 
plasma model and several gaussians is fitted to the 
baryon dominated non-blob reference region as shown in 
Figure [3] for Abell 520. No physical quantities are de- 
rived from this model, it is just designed to fit the data 
with a reduced x 2 ~ 1 for the numerical analysis of the 
spectrum (however, the resulting physical parameters are 
typical for clusters). All these baryonic model parame- 
ters, except the normalization, are frozen and the model 
describing the axion flux is added (the red line in figure 
[1]) . The continuous spectrum of decaying axions is repre- 
sented by redshifting the analytic expression in Equation 
[T]according to the distance of the cluster of galaxies. The 
composite model is fitted to the 2.0 — 9.5 keV spectrum of 
the dark matter blob region with the normalization of the 
two components as the only free parameters. There are 
no outstanding differences between the spectrum from 
the dark matter blob region and the baryon dominated 
reference region. 

An upper limit on the X-ray flux from the dark matter 
blob region is determined using two different methods. 
In the first method, all of the flux received from the re- 
gion is taken as a conservative upper limit on the flux 
originating from decaying dark matter. This method is 
independent of the assumed form of the decaying dark 
matter spectrum. The 3<r upper limit (with respect to 
the fitting residuals) on the flux is determined from the fit 
with the normalization given by the fitted value plus 3er. 
The obtained upper limits on the dark matter related 
X-ray luminosities are given in Table 1. In the second 
method, the decaying dark matter spectrum is assumed 
to have the form shown in Figure [T] (red line). An upper 
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FIG. 3: The observed spectrum of the blob region (red, 
squares) of Abell 520 and the spectrum of the reference region 
(black, triangles) with the fitted "basis model" (black, solid, 
reduced \ 2 ~ !)■ 



limit on the X-ray flux from a spectrum with this form 
which can be added to the "basis model" , is determined 
by increasing the derived model normalization systemat- 
ically in small steps while leaving the normalization of 
the "basis model" as a free parameter. When the re- 
duced x 2 reaches a value of 10 (A% 2 = 9), the 3er upper 
bound on the luminosity of the axion spectrum alone is 
determined. 



V. MASS DETERMINATIONS 

For Abell 520 the masses of the regions were deter- 
mined from weak gravitational lensing. The values are 
based on measuring the overdensity in the regions with 
respect to the mean density in a surrounding annulus 
with inner and outer radius of 0.85 arcmin and 4 arcmin, 
respectively. Hence, the mass value can be regarded as a 
conservative lower limit on the mass contained within the 
region. A detailed description of the data and methodol- 
ogy of the weak lensing analysis is given elsewhere [I?) • 
The mass of the Bullet Cluster blob region (SUB) is taken 
from 

The gas masses of the regions can be estimated as- 
suming a geometry for the gas distribution, e.g. that the 
gas is spherically distributed and with an average density 
derived from the thermal plasma model (MEKAL [ill ]) 
fitted to the spectra. As seen in Table 1, the gas masses 
of the blobs are much lower than their total masses (in 
agreement with the generally observed gas mass fraction 
of fgas ~ 0.11 [19]). Hence, the gas masses can be ne- 
glected and the mass of the dark matter is taken to be 
the total mass of the blobs determined from gravitational 
lensing. 
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Abell 520 




Bullet Cluster 






Value 


Ref. 


Value 


Ref 


Redshift 


0.2 


f21] 


0.29 


M 


D A , [Mpc] 


662 




872 




Region radius, 86 [arcmin] 


0.85 




0.66 




Chandra observation id 


4215 


[22] 


5356 and 5361 


[22] 


Exposure time [ksec] 


67 




179 




Blob region values 




M total [1O"M ] 


6.7 ±2.1 




5.8 


M 


M gas [10 13 M Q ] 


0.4 




0.7 


[14] 


Basis model, reduced \ 2 reference region 


0.75 




1.05 




Basis model, reduced x 2 blob region 


0.78 




1.04 




Basis + expected, reduced \ 2 


0.77 




1.05 




3<t total luminosity upper limit [10 44 erg/sec] 


0.2 




0.9 




3cr form dependent luminosity upper limit [10 44 erg/sec] 


0.2 




1.4 





TABLE I: The obtained values as described in Section ITVl and IVl for a cosmology with h = 0.71, Qm = 0.3, Qa = 0.7 [TTj ] . 



VI. LIFETIMES 

A simple estimate of the lifetime of the dark matter 
particles can be derived from the observed luminosities: 
L = TEjN, where T is the decay rate, E 1 is the photon 
energy, and N is the total number of particles. In the 
general case of non-relativistic radiative dark matter two- 
body decays, E 7 = m/2 for a dark matter candidate of 
mass m, N — X Mom /tti, where X is the mass fraction 
of the dark matter made up of the considered candidate. 
It is worth mentioning that we are considering the case 
where the many individual decay modes produce a wide 
KK-axion bump, which has the shape of the red (lower) 
curve in figure [T] The lifetime then becomes: 



1 IX Mom 
T= f = ^L— 



(2) 



For the two clusters the observational 3d upper limits 
on the luminosities lead to a lower limit on the mean 
lifetime of the order of r > 10 23 sec. The strongest con- 
straint comes from Abell 520 and gives an upper limit 
°f r ^ 6 x 10 24 sec assuming all of the dark matter to 
be made up of one single candidate with a radiative two- 
body decay. Of course the assumed axions produced in 
the stars (as predicted for the Sun), do not necessarily 
have to constitute all of the dark matter in the clus- 
ters. Still they would be confined by the gravitational 
potential of the cluster dark matter and an enhanced 
signal from decaying massive axions or the like would be 
expected from dark matter dense regions (such as the 
blobs). If this type of axions for instance contribute 
to 2% of the total dark matter mass, the lifetime con- 
straint relaxes to r > 10 23 sec (which corresponds to 
5077 < 3 x 1CP 15 GeV" 1 for a mean axion rest mass of 
5 keV). This first rough estimate of the lifetime is actu- 
ally not significantly different from that derived for solar 
KK-axions of r « 1.25 x 10 20 sec for a mean axion rest 
mass of 5 keV [l|. 

Moreover, we compare the cosmic densities of energy 



associated with dark matter and electromagnetic radia- 
tion [20| | since it will give a crude estimate of the radia- 
tive lifetime of dark matter particles on cosmic scale. We 
assume that the decay of a single kind of dark matter 
particles is responsible for the main part of the present- 
day electromagnetic radiation of the Universe. The total 
fraction of post-stellar radiation densities from resolved 
sources in radio-microwave, far infrared, optical, and X- 
7-rays is d to t ~ 2.4 x 10" 6 20]. With the age of the 
Universe taken to be Tuniverse 



4.5 x 10 17 sec 



111. this 



gives an order of magnitude estimate of the lifetime, in- 
dependent of particle type, of: 



Qdmtu- 



niverse 



4 x 10 22 sec , 



(3) 



which is similar to the life times derived above. It is en- 
couraging that studying KK-axion decay on scales from 
the Sun, to clusters of galaxies, to the size of the observ- 
able Universe results in life time constraints of similar 
order-of-magnitude. 



VII. SUMMARY 

Dark matter blobs in merging clusters of galaxies are 
excellent laboratories for constraining fundamental prop- 
erties of dark matter candidates, like their lifetime. One 
promising dark matter constituent might have properties 
similar to the generic KK-axion, presumably produced 
inside stars and proposed to explain the origin of the 
quiet solar X-ray spectrum which has remained elusive 
for several decades. Trapped by the deep gravitational 
potential of clusters of galaxies some of the KK-axions, 
or the like, eventually decay isotropically into two X-ray 
photons and thereby contribute to the diffuse intraclus- 
tcr X-ray emission. Hence, the X-ray emission from dark 
matter blobs can be taken as a conservative upper limit 
on the assumed total KK-axion luminosity, leading to a 
lower limit on the KK-axion lifetime of r > 10 23 sec. 
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This is several orders of magnitude larger than the life- 
time derived for solar KK-axions (10 20 sec). However, if 
KK-axions constitute less than a few percent of the dark 
matter mass, the two life time constraints are similar. 
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